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ABSTRACT 

We present millimeter and centimeter wave spectroscopic observations of the H n region NGC 6334A. 
We have mapped the source in several transitions of CO, CS, and NH 3 . The molecular emission shows a 
distinct flattened structure in the east-west direction. This structure is probably a thick molecular disk or 
torus (2.2 x 0.9 pc) responsible for the bipolarity of the near-infrared (NIR) and radio continuum emis- 
sion which extends in two “lobes” to the north and south of the shell-like H u region. The molecular 
disk is rotating from west to east (co & 2.4 km s" 1 pc” 1 ) about an axis approximately parallel to the 
radio and NIR emission lobes. By assuming virial equilibrium, we find that the molecular disk contains 
~2000 M©. Single-component gas excitation model calculations show that the molecular gas in the disk 
is warmer and denser (7^ ^ 60 K, n » 3000 cm" 3 * ) than the gas to the north and south (7^ ^ 50 K, 
n & 400 cm” 3 ). 

High resolution (~5") NH 3 (3, 3) images of NGC 6334A reveal several small (~0.1 pc) clumps, one of 
which lies southwest of the radio continuum shell, and is spatially coincident with a near-infrared source, 
IRS 20. A second NH 3 clump is coincident with an H 2 0 maser and the center of a molecular outflow. 
The dense gas tracers, CS J — 5 -» 4 and 7 -► 6, peak near IRS 20 and the H 2 0 maser, not at NGC 
6334A. IRS 20 has a substantial far-infrared (FIR) luminosity L ^ ~ 10 5 L 0 , which indicates the pres- 
ence of an O 7.5 star but has no detected radio continuum (F 6cm < 0.02 Jy). The combination of dense 
gas, a large FIR luminosity and a lack of radio continuum can best be explained if IRS 20 is a protostar. 
A third clump of NH 3 emission lies to the west of IRS 20 but is not associated with any other molecular 
or continuum features. The star formation activity in the region has moved west of NGC 6334A to IRS 
20 and the H 2 0 maser position. We suggest that NGC 6334A, IRS 20, and the H 2 0 maser spot are part 
of a “ protocluster ” of stars which is condensing from the massive molecular disk. The similarity between 
the structure around NGC 6334A and other large (r ~ 1 pc), massive (A/ ~ 10 3 M 0 ), rotating disks 
(K3-50A and G10.6 — 0.4) suggests that this may be a common mechanism by which open clusters form. 
Subject headings : ISM: H II regions — ISM: individual (NGC 6334A) — 

ISM: kinematics and dynamics — ISM: molecules — radio lines: ISM — 
stars: formation 


1. INTRODUCTION 

NGC 6334 is a southern H n region/star-forming cloud 
with at least seven sites of active massive star formation 
(d = 1.7 kpc; Neckel 1978). One of these sites of massive star 
formation, NGC 6334A, 2 has a shell-like morphology in 6 
cm radio continuum emission (Rodriguez, Canto, & Moran 
1982), with fainter radio emission extending in bipolar lobes 
to the north and south (Rodriguez, Canto, & Moran 1988, 
hereafter RCM88). NGC 6334A also has a bipolar mor- 
phology in near-infrared continuum emission (Harvey, 
Hyland, & Straw 1987), with extended nebulosity nearly 
coincident with the radio lobes. IRAS HIRES 12 and 25 /im 
data show double lobes to the north and south of NGC 
6334A as well (Fig. 1). Also associated with NGC 6334A are 
a CO “hot spot” (Dickel, Dickel, & Wilson 1977; Phillips, 
de Vries, & de Graauw 1986; this work), a Herbig-Haro- 
like object (Gyulbuaghian, Glushkov, & Denisyuk 1978; 
Bohigas 1992), high-velocity H 2 0 masers (Rodriguez et al. 
1980; Moran & Rodriguez 1980), and [C ii] 158 /xm, [O i] 
145 /x m, and [O i] 63 /xm emission (Kraemer et al. 1994). A 
far-infrared continuum source (X » 70 /xm, McBreen et al. 
1979; Loughran et al. 1986) coincides with the radio source. 

1 Departamento dc Astronomia, Universidad dc la Republics Monte- 
video, Uruguay. 

2 After the nomenclature of Rodriguez et al. (1982). NGC 6334A is 

spatially coincident with the FIR source NGC 6334IV (McBreen et al. 

1979). 


Two unresolved near-infrared (2.2 /im) sources (IRS 19 and 
IRS 20) are located between the extended continuum lobes 
(Harvey et al. 1987; Harvey & Gatley 1983, hereafter 
HG83). 

RCM88 suggested that the H ii region is confined in the 
east-west direction by a flattened structure of gas and dust, 
which would allow radiation to escape easily only to the 
north and south. There is some indirect evidence in support 
of this suggestion. A flattened dust structure oriented with 
its major axis in the east-west direction is implied by the 
elongated morphology seen in a map of 100 /xm optical 
depth (HG83). Moreover, De Pree et al. (1995) found an 
east-west velocity gradient in the H76a recombination line 
which they suggest arises from the ionized portion of a 
rotating accretion disk. 

We have imaged NGC 6334A in several millimeter and 
centimeter molecular lines in order to examine the distribu- 
tion and physical conditions of the molecular gas. The CO 
lines give the overall morphology of the gas, while the CS 
lines show where the dense gas lies. The NH 3 (3, 3) image 
reveals the small-scale structure of the warm dense molecu- 
lar gas near the continuum sources. We find that NGC 
63 34 A is surrounded by a massive, flattened, rotating struc- 
ture of molecular gas, just as RCM88 proposed. The emis- 
sion from the dense gas tracers does not peak at NGC 
6334A. Instead, they peak to the west, near the infrared 
source IRS 20 and an H 2 0 maser, which are probably pro- 
tostars. Apparently, the star formation activity has moved 
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Fig. 1. — Continuum emission at NGC 6334A: (a) 12 /mi IRAS HIRES image, peak — 1.4 kMJy sr“ contours at 200 MJy sr“ 1 intervals, lowest contour 
300 MJy sr - 1 . Positions of additional features are marked: 6 cm position (crass; Rodriguez et al. 1982), FIR peak (filled diamond ; McBreen et al. 1979), IRS 
20 (plus ; Harvey et al. 1987), and the H 2 0 maser (filled circle; Rodriguez et al. 1980; RCM88). (b) 25 pm IRAS HIRES image, peak = 9.5 kMJy sr -1 , 
contours at 800 MJy sr - 1 intervals, (c) 6 cm radio image (RCM88). (<t) K-band image (HG83). 


away from NGC 6334A. We suggest that a “protocluster” 
of stars is condensing out of the rotating cloud fragment. 

2. OBSERVATIONS 

We used the Caltech Submillimeter Observatory to map 
NGC 6334 in the CO J = 2 1 and 3 -► 2, the 13 CO 2 -► 1, 

and the CS 5->4 and 7 ->6 transitions. The CO 2 — ► 1, 
13 CO 2 — ► 1, and CS 5 -► 4 observations were made on 1993 
August 7-13. The CO 3 -+ 2 and CS 7 -+ 6 lines were 
observed simultaneously in dual-sideband mode, with one 
line in each sideband, on 1994 June 29-July 1. The facility 


230 and 345 GHz SIS receivers were used with the AOS 
backend (500 MHz bandwidth, 1024 channels). We imaged 
the cloud in CS 3 -> 2 at the NRAO 12 m telescope 3 at Kitt 
Peak on 1995 October 7. The A2 mm SIS receiver was used 
with the analog filterbank in the 256 channel, 1000 kHz 
channel -1 setting. All maps employed the “on-the-fly” 
mapping technique. At the CSO, individual grids (9 x 13 

3 The 12 m and VLA are operated by the National Radio Astronomy 
Observatory, which is operated by Associated Universities, Inc., under a 
cooperative agreement with the National Science Foundation. 
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and 13 x 19 pixels for 1993 and 1994, respectively) were 
mosaicked to create the larger maps. At the NRAO 12 m, a 
64 x 64 pixel grid was observed. A single spectrum of CO 

1 -► 0 emission toward NGC 6334A was obtained with the 
NRAO 12 m on 1995 April 5. Observational parameters are 
summarized in Table 1. 

The CSO data were processed with the CLASS and 
GreG packages and the NRAO data with the AIPS 
package. For calibration we used the standard chopper 
wheel method (Ulich & Haas 1976). Pointing was estab- 
lished through observations of Jupiter and varied by less 
than 5". The spectra were individually inspected and linear 
or parabolic baselines were removed as necessary. All data 
are presented on the T mb scale, that is, T mb = Tyrj mb (Table 
1). Intensities are estimated to be accurate to % 30%. 

The NH 3 (3, 3) observations were made with the VLA 3 , in 
the DnC array, on 1993 October 19 and 21. The NH 3 (3, 3) 
data were edited and calibrated in AIPS. The “pseudo- 
continuum ” data, from the “ channel 0 ” (the central 75% of 
the passband) were CLEANed and self-calibrated. A uni- 
formly weighted image was made in order to examine the 
continuum at high angular resolution. The continuum self- 
calibration solutions were then applied to the u , v line data. 
These data were used to produce a naturally weighted 
image cube, with a 30 kX taper applied. Emission-free chan- 
nels were averaged and subtracted from the data to produce 
a data cube which contained only line emission. 

3. THE MOLECULAR DISK 
3.1. Morphology 

The presence of a flattened structure of gas and dust was 
suggested by RCM88 to attenuate the radiation from the 
central star and to inhibit the expansion of the H n region in 
the east-west direction. Increased opacity in an east-west 
structure in the 100 /im optical depth map (HG83) provides 
indirect evidence in support of their suggestion. 

Our new observations show directly that the CO and CS 
emission is arranged in a flattened 2.2 x 0.9 pc structure 
(Fig. 2). In every transition, each of which is sensitive to 
different molecular gas densities, the emission is elongated 
in an east-west direction. The orientation of the major axis 
(P.A. » 90°) of the molecular gas emission varies by less 
than 10° among the various CO and CS lines. The FWHM 
axial ratio (Table 2) of the CO integrated intensity ranges 
from 1.6 for the CO 3 -► 2 emission to 2.5 for that of I3 CO 

2 -► 1. Along the minor (north-south) axis, the CS 3 -> 2 



Fig. 2. — Molecular emission at NGC 6334A. The (0, 0) position is 
a= 17 h 17 m 32 8 , <5— — 35°44'04", the position of NGC 6334 F. (a) CO 
2 — 1, peak integrated intensity = 708 Kkms' 1 , contours are at 50 K km 
s _l intervals. (fr) 13 CO 2 — 1, peak = 212Kkms~\ contours are at 15 K 
km s ~ 1 intervals, (c) CO 3 — 2, peak = 743 K km s ~ l , contours are at 50 K 
km s" 1 intervals. (7) CS 3 — 2, peak = 97 K km s“\ contours are at 
10 K km s“ 1 intervals; (c) CS 5 — 4, peak = 57 K km s~ contours are 
at 5 K km s _l intervals; (/) CS 7 — 6, peak = 88 K km s -1 , contours 
are at 10 K km intervals. Spectra along the horizontal line in (a) were 
fitted with Gaussian profiles for the velocity field diagram (Fig. 3) 

emission is barely resolved, and the CS 5 — ► 4 and 7 — ► 6 
emission are unresolved. Thus, the CS axial ratio of ~2 is a 
lower limit, and the dense gas, traced by the CS emission, 
may be more flattened than the more diffuse gas, traced by 
CO. A flattened structure similar to the molecular gas struc- 
ture we observe toward NGC 6334A arises naturally upon 


TABLE l 


Observational Parameters 


Species 

Transition 

Frequency 

(GHz) 

Velocity 
Resolution 
(km s -1 ) 

Beam 

(arcsec) 

Telescope* 


T 

<S 

CO 

7 = 1-0 

115.271 

2.60 

53 

KP 

0.55 

900 


J = 2-1 

230.538 

0.74 

26 

CSO 

0.74 

800 


7 = 3-2 

345.796 

0.49 

17 

CSO 

0.62 

4100 

13 co 

7 = 2-1 

220.399 

0.78 

27 

CSO 

0.74 

1750 

CS 

7=3-2 

146.969 

2.04 

42 

KP 

0.55 

650 


7 = 5-4 

244.936 

0.70 

24 

CSO 

0.74 

1000 


7 = 7-6 

342.883 

0.49 

17 

CSO 

0.74 

4100 

NH 3 

(7, K) = (3,3) 

23.870 

2.45 

6,6 x 5.1 b 

VLA 




1.3 cm continuum 

23.870 


2.3 x 2.0 C 

VLA 




» KP : NRAO 1 2 m at Kitt Peak ; CSO : Caltech Submillimeter Observatory ; VLA : Very Large Array. 
b Natural weighting with a 30 k2 taper. 
c Uniform weighting. 
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TABLE 2 

Properties of the Molecular Gas Structure 


Size 4 r h h Av M 

Line (arcsec x arcsec) (pc) (kms -1 ) (Af 0 ) 


CO 2 -► 1 126 x 54 1.1 3.9 2400 

l3 C0 2-l 96 x 38 0.8 2.5 730 

CO 3 ->2.. 39 x 24 0.3 2.4 280 

CS3-+2 86 x 46 0.7 2.6 710 

CS 5 -► 4 42 x 15 0.4 2.3 270 

CS7->6 34 x 17 0.3 1.5 90 


a Semimajor axis by semiminor axis at the half-power level of the 
integrated intensity. 

b At a distance of 1.7 kpc (Neckel 1978), 1' ~ 0.5 pc. 


the gravitational collapse of rotating clouds (cf. Boss 1987; 
Mestel 1965). 

3.2. Kinematics 

If the flattened structure of molecular gas arose from the 
collapse of a rotating cloud, we expect to see evidence of 
rotation in the kinematics of the gas. We fit Gaussian line 
profiles to the spectra along AS % — 471 //4 across the struc- 
ture (Fig. 2a), which are the spectra nearest to the major 
axis, and along Aa w —439", near the minor axis. Figure 3 a 
shows the central velocity from the Gaussian fits versus 
position across the molecular structure at A<5 « —471". In 
each transition, the velocity field shows a clear linear gra- 
dient of 2.4 km s" 1 pc - 1 along the major axis of the molec- 
ular emission. There is no systematic change along the 
minor axis from the average central velocity p lsr « — 3 km 
s" 1 . Figure 3 b shows the position-velocity diagram of the 

4 Offsets are with respect to the (0, 0) position at a = 17 h 17 m 32 s , 
<5 = — 35°44'04", the position of NGC 6334F. 


CO 2-^1 emission, which again shows the linear gradient. 
The gradient in the position-velocity diagram also shows 
that Gaussian profiles approximate the line shapes well. An 
east-west velocity gradient, that is, across the major (east- 
west) axis of the structure, is strong evidence that the molec- 
ular gas is rotating around a north-south axis. The com- 
bination of the linear velocity gradient, the elongation of the 
molecular emission, and the bipolarity of the diffuse contin- 
uum emission in a perpendicular direction leads us to con- 
clude that the structure is a rotating disk or torus of 
molecular gas and dust, with an angular velocity of co « 2.4 
km s' 1 pc -1 . This rotating molecular disk is the structure 
predicted by RCM88 which allows the continuum emission 
to escape only perpendicular to the structure. 

Unlike the velocity gradient seen on small scales (<30" 
from NGC 6334A) in hydrogen recombination lines (De 
Pree et al. 1995), in which the velocity increases to the 
northwest, the velocity of the molecular lines increases to 
the east. Central line velocities of the spatially integrated 
H76a and H92a profiles arer L sR~ —0.1 and— 0.32 kms -1 , 
respectively, and the velocity gradient across the H n region 
is ^75 km s" 1 pc -1 from southeast to northwest 
(P.A. « —60°; De Pree et al. 1995). De Pree et al. interpret 
this gradient as the combination of a north-south bipolar 
outflow and a rotating ionized disk with position angle 
P.A. « —90°. In contrast, the average central velocity of the 
molecular emission lies at p L sr ^ —3.0 km s -1 and has a 
position angle of P.A. % + 90°. Clearly the motions of the 
ionized gas are not coupled to the rotation of the molecular 
disk. This type of peculiar behavior, in which the ionized 
and molecular components seem to rotate in opposite direc- 
tions, has been reported in at least one other source, 
G34.3 + 0.2 (Garay & Rodriguez 1990), but there is no theo- 
retical explanation that we are aware of. We note that in 
NGC 6334A, the reddest H76a emission (De Pree et al. 
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Fig. 3. — (a) Velocity field diagram for the millimeter molecular lines. Velocities are taken from single Gaussian fits to each spectrum along a line of 
constant declination (A<5 « —471"; Fig. 1), approximately coincident with the major axis of the molecular disk. Note the smooth gradient across the disk. 
Error bars are formal errors to the Gaussian fits. (&) Position-velocity diagram of the CO 2 -* 1 emission across A<5 ss -471", the major axis of the disk 
(horizontal line in Fig. 2a). The coverage in right ascension is the same as in Fig. 3a. Contour levels are at 5% intervals of the 53 K peak (average) temperature; 
the lowest contour is at 50%. The higher velocity emission seen at a a 17M6 m 56 s is from the molecular outflow associated with the H z O maser (§ 5.2). 
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1995) appears coincident with NH 3 (3, 3) emission at the 
H n region (Fig. 4). This ionized gas may be the photoion- 
ized edge of the dense gas clump, not part of a rotating 
ionized disk. However, the apparent increase of the ionized 
gas velocity to the northwest, as opposed to the increase of 
the molecular gas velocity to the east, remains puzzling. 

3.3. Mass 

If the virial theorem holds, the mass of the disk can be 
estimated in each of the molecular transitions by 

rJ^Avf 
~ 0.4 G ’ 

where M is the mass and G is the gravitational constant 
(Binney & Tremaine 1987). The median radius, r h , which 
encloses half of the mass, was calculated by estimating the 
size of the major axis at the half-power contour of the inte- 
grated intensity maps. The numerical factor, 0.4, accounts 
for the difference between r g9 the radius enclosing all the 
gravitational mass, and r h , the median radius. The velocity 
change across the major axis, Ay, is determined by linear 
regression from the velocity field (Fig. 3a). 

Because of the high axial ratio in the dense gas tracers, we 
assume the disk is edge-on in the calculations of virial mass. 
If the disk is inclined by an inclination angle i, then the 
derived virial masses must by corrected by (sin i )~ 2 . By 
considering the case of a thin disk, we can place a limit on 
the inclination angle. If the disk were infinitely thin, then the 
observed thickness of the minor axis, r m{nor = 0.9 pc, would 
be due entirely to L The observed axial ratio sets a limit on 
i: sin i > (0.9/2.2) or 65° < i < 90°. For the lowest possible 
inclination, i v 65°, the virial mass increases by less than 
20%. Therefore, we conclude that inclination effects on the 
mass derivation are unimportant. 

The results of the mass calculations are given in Table 2. 
The virial masses derived vary from 90 M 0 for the CS 7 -► 6 
emission to 2400 M 0 for that of CO 2 -* 1. If the rotational 
velocities were due primarily to a central, condensed core of 
material, the calculated mass would be similar from each 


0 50 100 150 mJy 



RIGHT ASCENSION (B1950) 

Fig. 4. — NH 3 (3, 3) emission from NGC 6334A (contour interval is 
10% of the 0.82 Jy beam -1 kms' 1 peak), superposed on the 23 GHz 
continuum emission (peak flux = 150 mJy beam -1 ). Positions of addi- 
tional features are marked: 6 cm position (cross; Rodriguez et al. 1982), 
FIR peak ( filled diamond; McBreen et al. 1979), IRS 19 and 20 {plus ; 
Harvey et al. 1987), and the H 2 0 maser {filled circle; Rodriguez et al. 1980; 
RCM88). 


transition, not increasing. Also, the velocity field is linear 
across the disk, which is consistent with solid body rotation. 

The virial mass we derive from the CO 2 1 emission of 

M vir = 2.4 ± 1.1 x 10 3 M 0 is in excellent agreement with 
the molecular mass, M H2 = 2.2 + 1.4 x 10 3 M 0 , obtained 
by assuming a CO-to-H 2 conversion factor of 
N(H 2 )/I i-+q(CO) = 2.2 x 10 20 cm -2 (K km s' 1 ) -1 
(Combes 1991) and applying it to the mean CO (2 — ► 1) 
integrated intensity (averaged over a similar area). The 
uncertainty in the virial mass derivation is dominated by 
the uncertainty in Ay from the linear fit to the velocity field. 
For the molecular mass, the 30% uncertainty in the inten- 
sity calibration and an estimated 30% uncertainty in the 
CO-to-H 2 conversion factor contribute equally to the total 
uncertainty. Our use of CO 2 -> 1 instead of CO 1 -> 0 is 
unlikely to be problematic since the integrated intensity 
ratio in a 55" beam toward NGC 6334 A is CO (2 — ► l)/CO 
(1 -► 0) % 1, which indicates the CO is optically thick. These 
masses are slightly larger than that measured by Dickel et 
al. (1977), M H2 ~ 10 3 M 0 . Their measurements, however, 
were made at least 30" southeast of our CO peak, which 
probably accounts for the difference, 

To summarize, we find a large (2 pc diameter), massive 
(2000 M 0 ), rotating disk or torus of molecular gas sur- 
rounding NGC 6334A. The disk lies nearly edge-on in the 
east-west plane, is perpendicular to the lobes of radio and 
infrared continuum emission, and rotates about an axis 
parallel to the extended continuum emission. This rotating 
molecular disk is the structure of gas and dust which, based 
on the bipolarity of the continuum emission and the elon- 
gated 100 /xm opacity structure, RCM88 predicted was 
present. 

4. RELATION BETWEEN MOLECULAR GAS AND 

CONTINUUM EMISSION 

4.1. Physical Conditions 

In the CO 2 -► 1 line, the molecular gas apparentlyayoids 
the continuum lobes, as two distinct bubbles in the emission 
coincide with the radio continuum lobes (Fig. 5). RCM88 
suggested that the lobes of radio continuum radiation are 
caused by thermal emission from ionized gas and that the 
NIR emission is detectable because of lower extinction at 
the lobe positions relative to the surrounding region. We 
can test this hypothesis by determining the physical condi- 
tions of the molecular gas in the disk and toward the con- 
tinuum lobes. Qualitatively, comparison of the CO 2 1 

line emission with the 13 CO 2-> 1 emission (Fig. 6) shows 
that the CO opacity is indeed lower toward the lobes. 
Ratios of T mb (CO)/T mb ( i3 CO) indicate that t c0 « 40 in the 
disk, as opposed to r co « 15 toward the continuum lobes. 

A more detailed, multiline analysis confirms this qualit- 
ative result. We performed single-component model calcu- 
lations of non-LTE CO excitation to determine the physical 
conditions of the molecular gas. The model assumes that 
the emission originates in unresolved, homogeneous, spher- 
ical clumps. A photon escape probability function was 
included to account for the radiative excitation of optically 
thick lines (see Stutzki & Winnewisser 1985, and references 
therein). We modeled the emission from the first 1 1 levels of 
CO by varying the kinetic temperature, CO column density 
per velocity interval, and H 2 density until the best x 2 fit to 
the data was achieved. We used the collision rates of Flower 
& Launay (1985). We assumed that the beam-filling factor, 
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Fig. 5. — 6 cm radio continuum emission (contours ;RCM88) super- 
posed on the CO 2 -> 1 integrated intensity ( gray scale , 340-700 K km 
s~ *)■ Note the bubbles of decreased CO emission to the north and south of 
the peak near NGC 6334A coincide with the diffuse continuum emission. 



the fraction of the beam area covered by emission (<j> — 
T mb /7 Intrinsic)* was the same for each transition. A CO/ 13 CO 
abundance ratio of 60 was adopted. We did not include dust 
continuum emission, as it has negligible effect on low J CO 
excitation (cf. Jackson et al. 1995). Observed line parameters 
and the model results toward NGC 6334A, the northern 
lobe, located at (Aa, AS) « (— 410", — 415"), and the 
southern lobe, (Aa, AS) % ( — 417", —555"), are summarized 
in Table 3. 

By assuming a [CO]/[H 2 ] abundance ratio of 8 x 10“ 5 
(Frerking, Langer, & Wilson 1982) and multiplying the 
derived column density per velocity interval by the line 
width (Table 3), we find that the molecular hydrogen 
column density toward NGC 6334A is N H2 ~ 1.3 x 10 23 
cm -2 . The hydrogen column densities toward the northern 
and southern lobes are only N H2 ~ 4 x 10 22 cm -2 and 
N H2 ~ 9 x 10 22 cm -2 , respectively. This confirms our qual- 
itative estimate that the opacity in the disk is higher than 
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Fig. 6. — Map of optical depth as traced by the CO/ I3 CO 2 -* 1 inte- 
grated intensity ratio (gray scale), superposed on the CO 2 -► 1 integrated 
intensity map (contours). 


the opacity toward the continuum lobes. The lower column 
density toward the lobes supports the explanation of 
RCM88 that the diffuse NIR continuum radiation is detect- 
able because of lower opacity in the lobes as compared to 
that within the molecular disk. Further, the column density 
toward the southern lobe is roughly a factor of 2 higher 
than that toward the northern lobe, a result consistent with 
the suggestion of higher extinction toward the southern 
lobe (RCM88). 

If the expansion of the H n region is inhibited in the 
east-west direction by the molecular disk, the molecular gas 
in the disk must be denser than the surrounding gas. 
Further, when a cloud fragment collapses, as we propose 
the molecular disk has done, the density will increase as 
compared to density of the ambient parent cloud. Indeed, 
we find that the hydrogen volume density in the molecular 
disk is approximately 1 order of magnitude higher than that 
toward the lobes. This density enhancement is entirely con- 
sistent with both the confinement of the H n region and 


TABLE 3 


CO Line Parameters 


Parameters 

NGC 6334A 
(—416", — 471") a 

Northern Lobe 
( — 410" -415") 

Southern Lobe 
(-417", -555") 

T mb (K): 

CO 1 -► 0 

41 



CO 2 -► 1 

58 

36 

40 

CO 3 -► 2 

49 

26 

30 

I3 CO 2- 1 

27 

6 

14 

CO (2 - l)/CO (1 - 0) b 

1.4 



CO (3 - 2)/CO (2 -> 1) 

0.8 

6.7 

0.8 

( I2 CO 2- l)/( 13 CO 2- 1) 

3.7 

5.0 

4.7 

Model Results: 

TA K) 

60 

50 

50 

log n (cm 3 ) 

3. 3-3.8 

2.7 

2.5 

log (N co / Av) [cm “ 2 (km s - *) “ *] 

17.9-18.1 

17.7 

18.1 

log N H2 (10 22 cm -2 ) 

13 

4 

9 


3 The (0, 0) position is a = 17 h 17 m 32 s , S = -35° 44' 04", the position of NGC 6334F. 
b Integrated intensities were convolved to 55" (CO 1 -► 0 resolution) prior to the ratios. 
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p IG 7 fa) CO 2 -► 1 spectrum toward the H 2 0 maser position (—439, —457), with the Gaussian model fit and residuals. Residuals are on the same scale 

but are offset for clarity. Horizontal lines at 10 K show the velocity range over which the emission was integrated to produce Fig. lb. (b ) i Contour map of the 
integrated intensity of the CO 2 -+ 1 line wings. The redshifted emission ( dashed line) is integrated from 10-20 km s , and the blueshifted emission [solid line ) 
is integrated from -30 to -20 km s' 1 . Note that the outflow is centered on the position of the H 2 C> maser (filled circle), (c) Position-velocity diagram 
centered at the H 2 0 maser position (S = -35°51'40':4). Note the high-velocity emission at a % 17 h 16 m 56 B . 


with our scenario of a rotating molecular disk which con- 
densed out of the parent cloud. 

The CS emission was also modeled, but no valid solu- 
tions were found. That is, x 2 , the reduced x 2 , was never 
above the 5% confidence level. The CO model solutions, for 
comparison, had x? confidence levels of 55%-85%. We con- 
clude that the single-component model, while valid for the 
CO emission, does not adequately describe the CS emission. 


A two-component model, with a small, dense component 
embedded in an extended, diffuse component, might better 
describe the CS observations of NGC 6334A. 

In fact, large-scale mapping of the entire NGC 6334 
complex suggests that the CS 3 -► 2 emission and the CS 
5 4 and 7 -► 6 emission indeed trace distinct gas com- 

ponents (Kraemer & Jackson 1997). The CS 5-+4 and 
7 -> 6 emission appear only near the active star formation 
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sites. Apparently these lines specifically trace dense star- 
forming cores in NGC 6334. On the other hand, the CS 
3 -► 2 emission is more extended, similar to the CO 2 -► 1 
morphology. Thus, the CS 3 2 emission may be associ- 

ated with more diffuse gas than the gas traced by CS 5 -+ 4 
and J -> 6 emission. This is not unexpected, as the critical 
density of CS 3 -> 2 is 13 times smaller than that for CS 
7 -► 6. Unfortunately, a two-component model has more 
free parameters than we currently have observables. 
Further observations of CS isotopes or different CS tran- 
sitions are needed to properly constrain a two-component 
model. 

5. YOUNG STELLAR OBJECTS NEAR NGC 6334A 

5.1. IRS 20 

Two unresolved IR continuum sources, IRS 19 and 20 
(Harvey et al. 1987) (also known as IRS 2 and 3; HG83) lie 
between the NIR and radio continuum lobes (Fig. 1 c). IRS 
19 coincides with the southern rim of the radio shell, where 
there is a gap in the radio emission, and is probably associ- 
ated with the radio source. IRS 20, on the other hand, lies 
~20" west-southwest of IRS 19 (Fig. 4), just past the edge of 
the H n region. The two sources are strong at 20 //m (HG83) 
but not at J, H, or K bands. Indeed, IRS 20 was not even 
detected at J or H band ( J > 16.2, H > 14.0; Harvey et al. 
1987). 

A clump of NH 3 (3, 3) emission, which traces warm, dense 
gas, coincides with the position of IRS 20 (Fig. 4). The virial 
mass of this clump is M vir « 80 M 0 . Additionally, a number 
of other tracers of dense or photodissociated gas are associ- 
ated with IRS 20. For instance, the CS 5-* 4 and 7-»6 
emission, which are sensitive to densities of n> 10 6 cm -3 , 
peak not at NGC 6334A but instead at IRS 20 to within the 
positional errors (Figs. 2e and 2/). [C n] 158 jim emission, 
which traces photodissociated gas, is enhanced to the west 
of NGC 6334A, toward IRS 20, and the [O i] 145 fim 
emission peak is also coincident with IRS 20 (Kraemer et al. 
1994). The association of photodissociated gas with IRS 20 
is somewhat puzzling, as IRS 20 is not a radio source (Fig. 
4). In order to photodissociate CO and photoionize carbon 
without photoionizing hydrogen, 15 B1 zero-age main 
sequence (ZAMS) stars (Panagia 1973) are needed in a 
region ^0.04 pc across to be consistent with the observed 
luminosity (L^ %8 x 10 4 L 0 ; HG83). This required stellar 
density seems unreasonably high, and we suggest that the 
[C n] 158 jum and [O i] 145 pm emission are probably due 
to UV radiation escaping from the H n region. However, 
the infrared luminosity is far too high for IRS 20 to be 
merely a clump of gas illuminated by NGC 6334A. If we 
again assume a size of r « 0.02 pc, the size of the NH 3 
clump, and a projected distance of at least 0.2 pc from NGC 
6334A, IRS 20 intercepts less than 2% of the radiation from 
NGC 6334 A. The infrared luminosities of IRS 19, which is 
almost certainly the exciting source for NGC 6334A, and 
IRS 20 are both L, R - 10 5 L 0 (HG83). Clearly, IRS 20 must 
be self-luminous, and the CS and NH 3 trace dense circum- 
stellar gas surrounding IRS 20. 

IRS 20 has an infrared luminosity of L m « 8 x 10 4 L 0 , 
which requires a ZAMS star of O 7.5 (HG83). However, no 
radio emission has been detected from IRS 20 (F 6 cm < 0.02 
Jy ; RCM88), which puts a limit of B1 on the ZAMS spectral 
type. The combination of warm, dense gas, with high IR 
luminosity but no radio emission is exactly what is expected 


from a protostar. Thus, the evidence from our molecular gas 
observations supports the view (Harvey et al. 1987; 
RCM88) that IRS 20 is a site of active star formation near 
NGC 6334A. 

5.2. H 2 0 Maser 

Another NH 3 (3, 3) clump, M vir % 300 M 0 , coincides 
with an H 2 0 maser position (RCM88; Rodriguez et al. 
1980) approximately 15" north of IRS 20 and 20" west of the 
H n region (Fig. 4). The CO 2 -► 1 and 3-^2 emission both 
show extended line wings, indicative of a molecular outflow, 
toward the H 2 0 maser position (Fig. 7a). In the CO 2 -> 1 
emission, the red wings (10-20 km s _I ) and blue wings 
(—20 to —30 km s -1 ) separate into distinct lobes centered 
on the H 2 0 maser position (Fig. lb). The outflow can also 
be seen in the position-velocity diagram (Fig. 7c) at 
a « 17 h 16 m 56 s . We interpret these lobes as a molecular 
outflow centered on the H 2 0 maser. The gas in the outflow 
is at significantly bluer (t? LSR < — 20 km s -1 ) or redder 
(iYsr S: 10 km s* 1 ) velocities than the rotational velocities 
of the disk. Removing the contribution of the rotating disk 
from the postition-velocity diagram leaves the high-velocity 
gas at a % 17 h 16 m 56 s , with a much steeper velocity gradient. 
The combination of warm, dense gas with an H 2 0 maser 
and a molecular outflow, both signposts of star formation, 
indicates that this molecular clump may be yet another site 
of active star formation in the NGC 6334A area. The lack of 
IR emission, however, suggests that this source may be 
either less luminous or more deeply embedded than the IRS 
20 protostar. 

5.3. The Westernmost NH 3 Clump 

A third clump of NH 3 (3, 3) emission, M vir « 80 M 0 , 
~15" west of IRS 20, does not seem to be associated with 
any other molecular features. Also, there is no detected 
radio (RCM88; this work) or infrared (HG83; Harvey et al. 
1987) continuum emission at this location. Apparently, this 
clump of dense gas has not yet collapsed to form a proto- 
star. However, it does lie to the west of IRS 20, which is the 
direction in which star formation seems to be spreading 
from NGC 6334A. We suggest that although not currently 
active, this NH 3 clump may be the next site of star forma- 
tion within the molecular disk. 

6. CONDENSATION OF A PROTOCLUSTER? 

Molecular disks or toroids have long been thought to be 
the focusing mechanism of bipolar molecular outflows (e.g., 
Torrefies et al. 1983; Shu et al. 1993). However, when 
observed, these circumstellar disks generally have dimen- 
sions of at most a few tenths of parsecs, masses of M ~ 100 
M G or less, and are generally associated with low-mass star 
formation. A search of the literature finds three sources, 
K3-50A, G10.6 — 0.4, and V645 Cyg, with cloud properties 
similar to NGC 6334A; i.e., large (r ~ 1 pc), massive 
(M - 10 3 M g ), rotating disks. K3-50A, like NGC 6334A, is 
a shell-like H n region (Turner & Matthews 1984) with 
fainter radio emission extending to the northwest and 
southeast (De Pree et al. 1994, and references therein). 
K3-50A lies at the center of a large (d ~ 1 pc, M ~ 2 x 10 3 
M 0 ) molecular cloud (Vogel & Welch 1983) which rotates 
about an axis parallel to the extended radio emission. Ho, 
Terebey, & Turner (1994) found a massive (M ~ 10 3 M 0 ), 
rotating molecular core in G 10.6 — 0.4. This core is more 


622 


KRAEMER ET AL. 


Vol. 478 


compact (0.3 x 0.1 pc) and more flattened than the disks 
associated with NGC 6334A or K3-50A, but the mass esti- 
mate is similar. The case for a 10 3 Af 0 rotating disk in V645 
Cyg is less secure. There is some controversy concerning the 
correct distance to V645 Cyg, 3-6 kpc, which affects the 
estimates of the mass and size of the molecular gas struc- 
ture. The molecular outflow from V645 Cyg is oriented 
north-south at small scales (15"; Torrelles et al. 1987) but 
appears oriented northwest-southeast at scales of V and 
higher (Torrelles et al. 1987; Schulz et al. 1989). The 
ambient molecular gas at V645 Cyg is flattened east-west in 
the isotopic lines ( 13 CO 2 — ► 1, Torrelles et al. 1987; 13 CO 
and C 18 0 l->0, Verdes-Montenegro et al. 1991). The 
13 CO and C ls O l->0 position- velocity diagrams also 
show an east-west velocity gradient (Verdes-Montenegro et 
al. 1991). The elongation of the ambient gas and the velocity 
gradient perpendicular to the small scale molecular outflow 
are suggestive of a large (>1.1 x 0.5 pc, Verdes- 
Montenegro et al. 1991) rotating molecular disk. The mass 
estimate for the disk ranges from >140 M 0 (Verdes- 
Montenegro et al. 1991) to 2 x 10 3 M 0 (Torrelles et al. 
1987). 

Vogel & Welch (1983) first argued that the molecular disk 
of K3-50A might be a rotating protocluster. We suggest 
that these large, massive, rotating disks may, in fact, rep- 
resent a common phase in the formation of open clusters. 
These molecular disks are well matched to the size and mass 
of open clusters. Open clusters generally have masses of 
Afciuster ~ 250 M q (Binney & Tremaine 1987), which would 
require a typical star formation efficiency of ~ 10% for a 
2000 M q cloud. The molecular disks are smaller than the 
open clusters by approximately a factor of 2 (r h ~ 2 pc), but 
this may be due to dynamical relaxation of the system with 
time. 

The fact that a number of these large disks have now been 
discovered suggests the following mechanism for the forma- 
tion of open clusters. A slowly rotating molecular cloud 
fragment flattens because of its angular momentum (and 
possibly its magnetic field). An O star condenses in the 
center near the rotation axis, where angular momentum 
effects are minimized. An H ii region forms, but the stellar 
radiation and the expansion of the H n region are inhibited 
in the plane of rotation by the presence of the surrounding 
disk. The radiation from the massive young star can escape 
easily in the polar direction. This leads to the formation of 
two lobes of extended continuum emission perpendicular to 
the confining disk or torus of molecular gas and dust. If the 
stellar rotation axis is the same as that of the disk, and the 
jitellar outflow is directed poleward, the surrounding molec- 
ular gas and dust may not be immediately dispersed. Turb- 
ulence in the cloud fragment or shock waves from stellar 
winds may induce further star formation within the disk. 
This additional star formation, which propagates from the 
center outward, then leads to a gravitationally bound star 
cluster. 

The regions around NGC 6334A and K3-50A fit this 
scenario extremely well. In NGC 6334A, the H ii region has 
a shell-like morphology with fainter radio continuum emis- 
sion extending in bipolar lobes away from the shell. The H n 
region is su rr ounded by a large, massive disk of molecular 
gas which rotates about an axis parallel to the extended 
continuum lobes, NGC 6334A is at the kinematic center of 
the disk, which is where the first massive star is expected to 
form. The two protostar candidates, IRS 20 and the H 2 0 


maser, lie at the periphery of the radio shell, as if their 
formation was triggered by a shock from NGC 6334A. The 
third NH 3 clump lies further away from NGC 6334A and is 
apparently not (yet) an active star formation site. Thus, the 
star formation activity has moved outward from NGC 
6334A to IRS 20 and the H 2 0 maser but has not yet spread 
throughout the molecular disk. K3-50A also has faint, 
bipolar radio continuum lobes which extend from a shell- 
like H n region. Perpendicular to the lobes, a massive disk 
of molecular gas rotates around the H n region. The radio 
flux density is consistent with a late O star but is too low to 
account for the infrared luminosity L IR « 2 x 10 6 L 0 from 
the cloud. This suggests that the disk around K3-50A, like 
NGC 6334A, contains additional energy sources, that is, 
protostars. 

There remain a number of difficulties in this scenario of 
cluster formation. First, open clusters do not seem to rotate. 
This implies that the rotation of the disk must be slowed in 
some way. As Vogel & Welch (1983) suggest, when the gas 
from which the cluster formed disperses, the cluster will 
expand. As it expands, the rotational velocity will decrease. 
Additionally, if the open cluster forms in the inner part of 
the disk, angular momentum is minimal. Second, there is 
the question of star formation efficiency [SFE = 
M slars /(M stars + M gas )]. High-mass star-forming regions 
such as NGC 6334 tend to have efficiencies of - 5% or less. 
The SFE needed to form a gravitationally bound cluster is 
of order 50% (Lada & Lada 1991). However, if one con- 
siders only the gas involved in massive star formation (that 
is, the disk), and not the entire cloud (the NGC 6334 
complex) as is usually done, the SFE rises significantly. 
Further, it has been suggested (e.g., Bally & Lada 1991) that 
the birth of an O star marks the beginning of the end of star 
formation because of its disruptive influence on a cloud 
core. If the expansion of the H n region is inhibited by a 
molecular disk though, star formation may not be halted. 
These two effects may combine to raise the SFE within 
these molecular disks to an efficiency consistent with a 
bound cluster despite the presence of an O star. Finally, 
even if such protocluster disks are indeed common, they 
may be difficult to find. The distinctive edge-on geometry of 
the disks and bipolar continuum emission in NGC 6334A 
and K3-50A led to their discovery. If a cluster is forming in 
a disk which is face-on, the disk would not be easily recog- 
nized. A search for bipolar H ii regions might result in 
additional candidates. The molecular gas around these 
sources could then be examined for evidence of rotating 
disks. 

7. CONCLUSIONS 

1. We have found a massive (M ~ 2 x 10 3 M 0 ), flat- 
tened (2.2 x 0.9 pc) structure of molecular gas centered on 
the H ii region NGC 6334A. This disk or torus of gas and 
dust is rotating at an angular velocity of c u ~ 2.4 km s -1 
pc -1 about a north-south axis, which is parallel to the 
extended continuum emission. The molecular disk is 
responsible for confining the H ii region in the east-west 
direction, as suggested by RCM88. The radiation from the 
H n region can thus only escape perpendicular to the disk, 
to the north and south, and thus forms the radio and infra- 
red continuum lobes. 

2. We performed single-component model calculations 
of non-LTE CO excitation to determine the physical condi- 
tions of the molecular gas. Comparison of the CO bright- 
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ness temperatures and integrated intensity ratios with the 
excitation models show that the disk gas is warmer and 
denser (7^ « 60 K, n « 3000 cm’ 3 ) than the gas to the 
north and south (T k « 50 K, n & 400 cm’ 3 ). As RCM88 
predicted, the column density toward the continuum lobes 
is less than that toward the disk. 

3. NH 3 and CS emission, which trace dense molecular 
gas, peak at the infrared source IRS 20, west of the H n 
region NGC 6334A. IRS 20 has a high IR luminosity with 
no radio continuum and coincides with a warm, dense 
clump of molecular gas. A second clump of NH 3 coincides 
with the position of an H 2 0 maser and the center of a CO 
outflow. The most likely explanation for these phenomena 
is that IRS 20 and the H 2 0 maser spot are protostars. A 
third NH 3 clump, west of IRS 20, has not yet formed a 
young stellar object but is likely to be the next site of star 
formation activity in the region. 

4. The molecular disk centered on NGC 6334A joins 
K3-50A, G10.6-0.4, and V645 Cyg as the fourth example 


of a large, massive, rotating disk. The similarities in molecu- 
lar and continuum properties between NGC 6334A and 
K3-50A (i.e., the massive rotating molecular disk and the 
perpendicular bipolar continuum emission) are remarkable. 
We suggest that such these massive, rotating disks provide 
natural environments for the formation of open clusters. 
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